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Abstract. The effects of rapid rotation and bi-stability upon the density contrast between the equatorial and polar directions 
of a B[e] supergiant are re-investigated. Based upon a new slow solution for different high rotational radiation driven winds 
(Cure 2004) and the fact that bi-stability allows a change in the line-force parameters (a, k, and 6), the equatorial densities are 
about lO'-lO"* times higher than the polar ones. These values are in qualitative agreement with the observations. 
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1. Introduction 

B[e] supergiants belong to a post-main sequence evolutionary 
stage of massive stars, their characteristics are; Balmer lines 
in emission, sometimes with P Cygni profiles with equivalent 
widths for Ha greater than 100 A, low-excitation permitted 
emission lines, predominantly of singly ionized metals, forbid¬ 
den emission lines of [O I] and [Fe II] and strong near/mid 
infrared excess due to hot circumstellar dust, indicating dust 
temperatures of lOOOK (Zickgraf et al. 1986, 1992). 

A large percentage of B[e] supergiants (~ 70%-80%) show 
a hybrid spectrum (Zickgraf et al. 1985,1986) that comes from 
the simultaneous observations of optical narrow low-excitation 
emission lines (~ 100 kms^^) and broad UV high-excitation 
absorption lines (C IV, Si IV and N V) with terminal velocities 
similar to early B supergiants winds, of order of ~ 1000 km . 
An empirical model that describes these hybrid spectral char¬ 
acteristics was suggested by Zickgraf et al. (1985) in terms of 
a two-component stellar wind, consisting of a fast radiation- 
driven wind (Castor et al. 1975, hereafter CAK) from polar 
latitudes and a slow and dense expanding disk in the equato¬ 
rial regions of the star. This two-component wind has been 
confirmed from medium resolution spectropolarimetry for the 
evolved B[e] star HD 87643 (Oudmaijer et al. 1998). 

The bi-stability mechanism induced by rotating radiation 
driven winds was introduced by Tamers & Pauldrach (1991) 
to explain the formation of outflowing disks around early-type 
stars. This mechanism, which causes a drastic change in the 


wind structure, was thought to be related to the behaviour of 
the Lyman continuum optical depth, ti, with the stellar lati¬ 
tude. For tl < 1 the wind is fast but for t/, > 3 the wind is 
slow. The bi-stability manifests itself at certain spectral type 
or temperature. Tamers et al. (1995) have determined terminal 
velocities, Voo, of a large sample of early-type stars concluding 
that the bi-stability jump is present around Tefr = 21,000 K. 
Vink et al. (1999) have theoretically shown that the bi-stability 
mechanism is located at Teff = 25,000 K and that it is mainly 
due to the radiative acceleration by iron, caused by the recom¬ 
bination of Fe IV to Fe III. 

Pelupessy et al. (2000) have calculated density contrast (ra¬ 
tio between equatorial and polar densities) in a B[e] supergiant 
for rotationally induced bi-stability models applying multi¬ 
scattering line-force parameters above and below the critical 
temperature of the bi-stability jump. They show that the ra¬ 
tio between equatorial and polar densities is about ~ 10, and 
state that this value is a factor 10 times smaller than Bjorkman’s 
(1998) calculations. Porter (2003) also asserts the difficulties in 
modelling optical-near-IR emission for B[e] supergiants with 
disk model’s density structure produced by either a bi-stability 
wind or a Keplerian viscous disk. 

Pelupessy et al. (2000) wind solutions were computed con¬ 
sidering values of Q = Vrotjvhrkup ^ 0.6, where Vrot is the 
equatorial rotational speed and Vhrkup is the break-up speed. 
However, B[e] supergiants are located near the Eddington limit 
(Zickgraf et al. 1986). Here, critical rotation speed is reached 
at a much lower star’s rotational speed. Velocities of about 
200 km would make the star rotate sufficiently close to the 
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break-up speed to produce observable effects. Consequently, 
the wind characteristics near the equator are expected to dif¬ 
fer from the polar wind. Langer (1998) proposes the fl- 
mechanism for these type of stars. His model calculations sug¬ 
gest that despite the loss of angular momentum due to mass loss 
and increasing radius a rapidly rotating massive main-sequence 
star could remain for a substantial fraction of its lifetime close 
to the so-called Q-limit which designates the limit of critical 
(or break-up) rotation. 

In order to investigate the influence of rotation in radiation 
driven winds (even for a star rotating up to break-up rotational 
speed), Cure & Rial (2004) performed a topological analysis of 
the rotating CAK model, finding that the line-force parameter 
d (Abbott 1982), that accounts for changes in the ionization 
of the wind, leads to a bifurcation in the solution topology 
as a function of the rotational speed, shifting downstream the 
location of the critical point of the CAK (x-type) singular 
point. Thus a higher mass-loss rate and lower terminal velocity 
wind is attained in comparison to the frozen-in ionization 
(6 - 0) case. 

Furthermore, Cure (2004) proves that the standard solution 
(hereafter the fast solution) of the m-CAK wind model (Friend 
& Abbott 1986, Pauldrach et al. 1986) vanishes for rotational 
speeds of ~0.7 - 0.8 Vbrkup, and there exists a new solution, that 
is much denser and slower than the known standard m-CAK 
solution. We will call it hereafter the slow solution. 

The purpose of this work is to re-investigate the forma¬ 
tion of an equatorial disk-wind for rapidly rotating B[e] su¬ 
pergiants, taking into account; 1) the fast and slow solutions of 
rotating radiative driven winds that depend on the assumed ro¬ 
tational speed Q and 2) bi-stability line-force parameters. 

In section 12 we discuss the adopted line-force parameters. In 
section 13 we show results for CAK and m-CAK models with 
fast and slow solutions. Section |4] is devoted to the analysis 
of the influence of changes in ionization throughout the wind. 
Discussion and conclusions are presented in sections |3and|6l 
respectively. 

2. Rotating bi-stability parameters 

In order to investigate the influence of the rotation and the bi¬ 
stability jump in forming a disk-wind, we solve the non-linear 
momentum equation for the CAK wind and for the m-CAK 
wind, in both polar and equatorial directions. Details and cal¬ 
culation methods here used for the CAK wind are found in Cure 
& Rial (2004) and for m-CAK in Cure (2004). 

Prior to solving CAK and m-CAK momentum equations, 
we have to know the line-force parameters. The parametriza- 
tion from Abbott (1982) of the line-force uses 3 parameters (a, 
k and 6). These parameters are obtained from the fitting of the 
net line acceleration of hundreds of thousands NLTE line tran¬ 
sitions computed by solving consistently the radiative transfer 
and hydrodynamics. A different approach for calculating the 
line-force parameters that include multi-line effects for tem¬ 
peratures above and below the bi-stability jump was developed 
by Vink et al. (1999). They use the Monte-Carlo technique and 
assume a /3-field for the velocity law. 


Table 1. Bi-stability line force parameters 


T[K] 

a 

k 

5 

30,000 

0.65 

0.06 

0 

17,500 

0.45 

0.57 

0 


We want to stress, that no calculation of the line-force pa¬ 
rameters has been performed for the slow solution and this is 
beyond the scope of this study. Therefore, in order to solve the 
wind momentum equation, we adopt Pelupessy et al. (2000) 
line-force parameters, a and k, since they have been calculated 
for both sides of the bi-stability jump. These line-force param¬ 
eters are summarized in Table[T] 

3. Wind model results 

In order to compare our results with the ones from Pelupessy 
et al. (2000), we adopt the same B[e] supergiant star they used; 
Tgg — 25,000 K, MIMq - 17.5, L/Lq - 10^ and solar abun¬ 
dance. For the lower boundary conditions for polar and equato¬ 
rial directions, we use the same procedure as Stee & de Araujo 
(1994), i.e., after solving the momentum equation (CAK and 
m-CAK) in the polar direction, with the surface boundary con¬ 
dition, T = 2/3 (electron scattering optical depth), we obtain 
the value of the polar photospheric density, ppiRt). This value 
of the photospheric density is then used as surface boundary 
condition for the equatorial direction, for both CAK and m- 
CAK wind models, respectively. 

In this work we have not taken into account neither the 
change in the shape of the star nor gravity darkening (von 
Zeipel effect) nor the modification of the finite-disk correction 
factor due to the rotation (Cranmer & Owocki 1995 equation 
[26], Pelupessy et al. 2000). However, we expect that these ef¬ 
fects may hava a small influence on the fast solution (see sec¬ 
tion |^^). The study of their influences in the slow solution 
will be the scope of a forthcoming article. 

3.1. The rotating bistability CAK wind 

As a first step to study the combined effect of bi-stability and 
high rotation speed we investigate the CAK wind. We solve the 
isothermal (Teff = 25,000) CAK wind for our B[e] supergiant, 
for the polar(equatorial) direction we use the line-force param¬ 
eters determined for temperatures above(below) the bi-stability 
jump, see Table [2 

Figure [2 shows the density profile in both directions, po¬ 
lar (dotted line) and equatorial (continuous lines), for different 
values of the Q parameter and Figure |2 shows ratios between 
the equatorial, p^, and polar, pp, densities. Our result, without 
finite-disk correction factor, gives a density contrast similar 
to the one obtained for Q = 0.6 by Pelupessy et al. (2000), 
who use the finite disk correction factor. Density contrast for 
Q > 0.7 are a few times denser. While the CAK model can 
be calculated for any rotational speed, Pelupessy et al. (2000) 
could not calculate m-CAK models for D. > 0.6, the reason 
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Table 2. Parameters of the calculated CAK models: ter¬ 
minal velocity, Voo (kms^^), the local mass loss rate 
(10^^ MQyr^^) and the location of the critical point. 



Q 

Fee 

F„ 

rdR. 

pole 

0.0 

826 

0.181 

1.60 

equator 

0.6 

379 

0.718 

9.83 

equator 

0.7 

348 

0.730 

11.47 

equator 

0.8 

319 

0.743 

13.13 

equator 

0.9 

293 

0.756 

14.80 

equator 

0.99 

273 

0.770 

16.33 



Fig. 1. CAK model: density (in grcm versus r/R, - 1. 
Polar density is in dashed-line; equatorial densities are in 
continuous-line, the higher is Q., the higher is the density 
(Q = 0.6,0.7,0.8,0.9,0.99). 



Fig. 2. CAK model: density ratio pe(Q.)lpp versus r/R, - 1. 
The higher is Q, the higher is the density contrast. Curves are 
for Q = 0.6, 0.7, 0.8, 0.9, 0.99, respectively. 

is because the fast solution ceases to exist (Cure 2004). From 
Figure 12] we can see that a disk wind structure is formed, the 
higher is Q, the denser and larger is the disk close to the photo¬ 
sphere and up to r ~ 2/?,, region where the equatorial density 
reaches (as function of Q) values of hundreds times the value 
of the polar density. Far from the stellar surface and up to hun¬ 
dreds stellar radii, the ratio pdPp reaches values of ~ 10. 
Hence, a rotating bi-stability radiation driven wind forms a 
disk at the equatorial latitudes of these stars. Table |2| summa¬ 
rizes the polar and equatorial values of the terminal velocity. 



Fig. 3. m-CAK model: density (in grcm ^) versus r/f?, - 1. 
Polar density is in dotted-line; equatorial density for Q = 0.6 
(fast solution) is in dashed-line and equatorial densities for Q = 
0.7, 0.8, 0.9, 0.99 are in continuous-line, the higher is Q, the 
higher is the density. 



Fig. 4. m-CAK model: density contrast versus r/f?, - 1, 
dashed-line is for Q. - 0.6 and continuous-line are for Q = 0.7, 
0.8,0.9,0.99. The higher is Q, the higher is the density contrast. 


Fm, the ’local mass loss rate’ (see Pelupessy et al. 2000, equa¬ 
tions [10] and [11]) - the total mass loss rate if the solution 
for this latitude were valid for a spherical star - and the lo¬ 
cation of the critical (singular) point. The terminal velocities 
are higher than the observed ones in B[e] supergiants disks 
(Zickgraf, 1998). This fact might be due to the assumed line- 
force parameters. There are no calculations of these parameters 
for the CAK model with rotational speeds. Therefore, these re¬ 
sults represent our first approximation to modeling the outflow¬ 
ing disks of these objects. 

3.2. The rotating bistability m-CAK wind 

We present here the results of m-CAK wind models (see so¬ 
lution scheme in Cure 2004), taking the same set of line-force 
parameters given in Tableland different values of the rota¬ 
tional speed for the equatorial direction. Due to the fact, that the 
existence of the fast or slow solution depends on the rotational 
speed, we analyse the fast solution for Q. - 0.6 (~ upper limit 
of Q for fast solutions) and then for higher rotational speeds we 
obtain slow solutions. We show that a disk-wind with a large 
density contrast is formed for high rotational speeds. Our re- 
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suits are summarized in Table 0 displayed in figures 0 and 0] 
and will be discussed below. 

3.2.1. The fast solution 

In order to compare our results with the ones from Pelupessy 
et al. (2000), we solve the m-CAK momentum equation, for 
Q = 0.6. Note that, there are no slow solutions that satisfy the 
lower boundary condition for this value of Q. 

This fast solution is shown in figures 0 and 0] by dashed- 
line. Its density is lower than the densities from slow solutions 
(continuous-lines, see section l3.2.2> and higher than the polar 
density (Q = 0). We obtain a density contrast of about 10 for 
almost all the wind, similiar to Pelupessy et al. (2000) result. 
Since Pelupessy et al. (2000) included the effects of: change in 
the shape of the star as a function of the rotational speed, de¬ 
pendence of the temperature on the latitude and the finite disk 
correction factor due to an oblate star, we expect that these ef¬ 
fects have a small influence on the fast solutions. 

The increase in the density contrast in the region close to the 
photosphere (see figure© is due to the centrifugal force and the 
consequently higher mass-loss rate of the fast solution when 
rotation is included (Friend & Abbott 1986). 

3.2.2. The slow solution 

We have calculated slow solutions from the m-CAK momen¬ 
tum equations in the equatorial direction for rotational speeds 
for Q = 0.7, 0.8, 0.9, 0.99. Density profiles from figures^and 
© for both CAK and m-CAK solutions respectively, show a 
similar behaviour. The differences near the photosphere and up 
to ~ 2Rt are due to the finite disk correction factor, fo, but for 
radii larger than two stellar radii, both density structures are al¬ 
most the same. Near the photosphere fo is less than one, thus 
a lower mass loss rate is attained. Further out in the wind the 
value of fo is greater than one, so the plasma is accelerated to 
higher terminal velocities. 

Figure©shows the density contrast profile, which are larger 
than the corresponding CAK models (see figure |2ll. This is 
due to the inclusion of fo in the momentum equation, giv¬ 
ing a polar flow (Q - 0) less dense than the CAK case. 
Density contrasts reach values around thousand for radii less 
than ~ 2Rt and a value of hundred is maintained by the wind 
up to hundreds of stellar radii, almost independently of Q. This 
result concerning the disk behaviour is in qualitative agreement 
with the values estimated from observations by Zickgraf et 
al. (1985,1986,1992), Zickgraf (1998), Oudmaijeretal. (1998) 
and Bjorkman (1998). 

4. Changes in ionization throughout the wind 

Prinja et al. (2005) investigate the ionization structure of early- 
B supergiant winds and demonstrate that the wind ionization 
increases with distance from the star. This structure is different 
that the one exhibits by an O-star wind. Therefore in order to 
study the effect of changes in ionization with radial distance in 
a fast rotating supergiant, we explore the influence of the pa¬ 
rameter 6 in the m-CAK model. 


Table 3. Parameters of the calculated m-CAK models: ter¬ 
minal velocity, Vca (kms^^), F,„, the local mass loss rate 
(10^^ MQyr~^) and the location of the critical point. 



a 

Poo 

F,n 

rdR, 

pole 

0.0 

2287 

0.109 

1.027 

equator^' 

0.6 

lAl 

0.724 

1.057 

equator 

0.7 

340 

0.838 

12.75 

equator 

0.8 

312 

0.851 

14.21 

equator 

0.9 

289 

0.864 

15.07 

equator 

0.99 

266 

0.875 

17.14 


“ Fast solution values. 



Fig. 5. Density contrast for the m-CAK wind with 6^0. 
Parameters a and k are given in Table© 

Table 4. Parameters of the calculated m-CAK models with 6 + 
0 from Abbott (1982), terminal velocity, Voo {km i“'), the local 
mass loss rate, F,„ (lO^^MQjr^') and rdRt, the location of 
the critical point. 



Q. 

Too 

F„ 

rdR. 

pole 

0.0 

1475 

0.097 

1.053 

equator 

0.7 

318 

1.042 

15.07 

equator 

0.8 

292 

1.087 

16.04 

equator 

0.9 

267 

1.138 

18.41 

equator 

0.99 

248 

1.188 

19.88 


We adopt 6 values from Abbott (1982), for the polar direc¬ 
tion, (5 = 0.12 (T - 30,000) and for the equatorial direction, 
6 = 0.089 (T - 20,000). The a and k values are the ones from 
Table© Figure©shows density contrast profiles for Q > 0.7, 
which exhibit the same behaviour as in the frozen-in ioniza¬ 
tion case (section l5.2.2t and values of the same order. We have 
also calculated m-CAK models with the 6 line-force parame¬ 
ter from Shimada et al. (1994), arriving at the same conclusion, 
that the b-parameter has a small influence in the density con¬ 
trast value. Model parameters are summarized in Table© 
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5. Discussion 

We want to stress the importance of the combined effect from 
slow and fast solutions with bi-stability line-force parame¬ 
ters in forming an outflowing disk wind in B[e] supergiants. 
Density contrast of order of 10^ up to large distances from the 
star are attained. 

This theoretical value is in qualitative agreement with the val¬ 
ues derived from observations of order of 10^-10^ (Zickgraf 
et al. 1989, Zickgraf 1998 and references therein, Bjorkman 
1998). 

Previous simulations of disk formation in rotating radia¬ 
tion driven winds induced by bi-stability have underestimated 
the density contrast, mainly due to: a) The use of a /3-field 
(with = 1, Tamers & Pauldrach 1991) to describe the wind 
velocity profile, even for high rotational speeds where this 
approximation fails (Cure 2004). Table [1] from Tamers & 
Pauldrach (1991), shows that t/, > 3 exactly when the standard 
fast solution do not further exist, and b) Pelupessy et al. (2000) 
calculations, based in the fast solutions, were restricted to val¬ 
ues of D < 0.6 and these rotation values are not high enough 
for developing a dense disk. 

A dense disk is formed when the slow solution starts to 
exist. For our test star, this occurs for Q > 0.7. This condi¬ 
tion is in agreement with the estimation of 0.74 < Q < 0.79 
done by Zickgraf (1998) in order to reproduce observable ef¬ 
fects in the structure of stellar winds. However, observational 
rotation speeds of B[e] supergiants have high uncentainties, be¬ 
cause only a few stars show photospheric absorption lines ap- 
propiated for the measurements of V sin(/). The inferred obser¬ 
vational value of Q lies in the range 0.4-0.7 (Zickgraf 1998). 

We have verified that the <5 line-force parameter, which is 
related to changes in ionization troughout the wind, modifies 
the wind structure, mass-loss rate, terminal velocity and loca¬ 
tion of the singular point (compare Tables 0 and 0}. Despite 
these changes, the density ratio pdpp is almost the same as in 
the frozen-in ionization case. 

Furthermore our result for the density contrast based on the 
m-CAK fast solution for Q. - 0.6 are of the same order as 
Pelupessy et al. (2000), eventhough they included a more de¬ 
tailed description of the distorsion of the star due to the rota¬ 
tional speed. 

The existence of a fast solution for the polar direction and a 
slow solution for the equatorial direction, calculated using bi¬ 
stability line-force parameters, give a density contrast of order 
of 100. 

Since most of the B[e] supergiants in the H-R Diagram are 
located below the bi-stability jump temperature (25,000K), in 
our conception, the theoretical explanation for the existence of 
a two-component wind model (Zickgraf et al. 1985) is due to 
the nature of the solutions of rapidly rotating radiation driven 
wind. The change (jump) from fast solution to slow solution at 
some latitude accounts by itself a two-component wind, where 
each solution structure has its own set of line-force parameters. 
This picture would be remarked for cases when the bi-stability 
jump is present. 

Another important aspect to remark is the scarcity of self- 
consistent calculations of line-force parameters k, a, 6 for the 


m-CAK fast solution and the lack of calculations for our slow 
solution. The uncertainty in the values of the parameters re¬ 
flects in the value of the terminal velocity, mass loss rate, as 
well as in the density contrast. Specifically, the predicted ter¬ 
minal velocities, see tables|3and|4] are about 2-3 times greater 
than values inferred by observations (Zickgraf 1998). 

Therefore our results, that combine fast and slow wind solu¬ 
tions, are a first approximation to re-investigate disk formation 
in high rotating stars with radiation driven winds. A detailed 
wind model needs a self-consistent line-force parameter cal¬ 
culations for both, fast and slow solutions. 

6. Conclusions 

We have revisited radiative driven wind models for a high ro¬ 
tating B[e] supergiant (Q > 0.6) assuming a change in the line- 
force parameters due to the bi-stability jump. The existence of 
slow and fast solutions in the model, predicts density contrast 
which are of order of 10^ - near the stellar surface (r < 2 
R,), while outside it falls to values of about tens to hundreds 
and the disk extends up to ~ 100 stellar radii. Comparing the 
density contrast predicted by the CAK and m-CAK models, 
we conclude that the m-CAK model better describes disk for¬ 
mation in B [e] supergiants being in qualitative agreement with 
observations. 
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